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Abstract 

We study the ferromagnetic Ising model on the Sierpihski gasket (SG), where the 
spin-spin interactions depends on the direction. Using the renormalization group 
method, we show that the ratios of the correlation lengths restore the isotropy of 
the lattice as the temperature approaches zero. This restoration is cither partial or 
perfect, depending on the interactions. In case of partial restoration, we also evaluate 
the leading-order singular behavior of the correlation lengths. 

1 Introduction 

Restoration of microscopic isotropy in fractals with microscopic anisotropy is observed in 
various physical systems, including percolation [I], diffusion and random walks [2, 3J, resister 
networks [1], and Gaussian field theories [5]. Such phenomena are unique in the sense that 
they are absent in uniform media, while they are universal in the sense that they are observed 
in a wide class of fractals. 

The isotropic spin 1/2 Ising model on the Sierpihski gasket (SG) is known to have no phase 
transitions at finite temperature. But the correlation length exhibits a drastic divergence as 
the temperature approaches zero [3, [8]. It would be quite interesting to determine whether 
the anisotropic version of the model exhibits isotropy restoration. 

Indeed, Brody and Ritz [9j, and later Bhattacharyya and Chakrabarti pi)] studied this 
problem in a partially anisotropic Ising model (which is obtained by setting J 2 = J3 or 
Ji = J 2 in our model) on the SG. Interestingly, they reached opposite conclusions: Ref. [9] 
concluded that there is a complete restoration, while Ref. [10] concluded that restoration is 
impossible. 

In the present paper, we carefully examine this (controversial) problem of isotropy restora- 
tion in the Ising model on the SG. We study a general anisotropic model in which spin-spin 
interactions take three different values, according to the directions of the bonds. 

In order to test for the possible restoration of isotropy, we define the correlation length 
for each of the three directions of the lattice. When the temperature is sufficiently high, the 
ratios among the three correlation lengths are almost the same as those of the interactions. 
We conclude that the behavior of the correlation lengths respects the anisotropy of the 
system. 
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A nontrivial question is whether the isotropy is restored when the temperature approaches 
zero and the correlation lengths diverge. We have examined the behavior of the correlation 
lengths using the renormalization group, and found that the isotropy is indeed restored. 
However, we also found that, depending on the interactions, the restoration can be either 
perfect or partial. When perfect restoration takes place, all the ratios among the correlation 
lengths approach unity as the temperature approaches zero. When a partial restoration takes 
place, on the other hand, the ratio of the two correlation lengths approaches unity, and the 
other correlation length becomes infinitely larger, as the temperature approaches zero. 

It should be stressed that our result agrees with neither the conclusion of Ref. [9] nor 
that of Ref. [ID] . We have shown that the actual phenomena are much richer than those 
predicted in Refs. [9] or [TUJ . We do, however, understand why such different conclusions 
have been reached. As we discuss in Section 13.51 both Refs. [9] and [10] essentially rely on 
simple approximate forms of the RG equation. But, as will be clear from the careful study 
presented here, the RG flow in this problem cannot be reduced to a single approximate map; 
one needs to use different forms in different regions of the parameter space. We thus believe 
that our analysis, which is clearly more elaborate and careful than the previous two, resolves 
the controversy and reveals the true behavior of the model. 

We also stress that no distinction between perfect and partial restoration of isotropy, as 
we have found in the present model, has been observed in other models [H El El HI [5] in 
which only perfect restoration was found. We still do not know why partial restoration is 
found only in the Ising model. 

Let us stress that the restoration of isotropy is a phenomenon peculiar to physical systems 
on fractals. In the Ising model on regular lattices, it is believed that the ratios of the 
correlation lengths respect the anisotropy of the underlying interactions, no matter how 
large the correlation lengths are. 

2 Model and main results 

Here we define the model and state our main results. 

We construct the SG following the standard procedure (see Fig. [T]). We consider an 
equilateral triangle formed by bonds of unit length, and call this the first generation SG. We 
then put three of these triangles together to construct the second generation SG, which is 
identical to the original triangle with all sides doubled in length. In the next step, we put 
three second generation SGs together and form the third generation SG. By continuing this 
procedure recursively, we construct the n-th generation SG. 

We define an Ising model on the n-th generation SG by associating with each site an 
Ising spin. The Hamiltonian of this system is 

Tin = -^l(o"0QO"l0 + O"l0O"20 + o"oiO"n H ) 
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Figure 1: The first three generations of the Sierpihski gasket. 
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Figure 2: Ising model on the Sierpihski gasket. The spins are labeled according to the oblique 
coordinates. 

-J 2 ((Tiocroi + cr 2 ocrii + 0"nO"o2 H ) 

-^3(o"ooO"oi + cioCu + o"oiO"02 + •••)> (!) 

where the sums are taken over all the pairs of neighboring spins. We label each spin variable 
according to the oblique coordinates, as in Fig. [2j The positive constants Ji, J 2 and J3 
denote the exchange interactions in the three directions. 

The thermal expectation value of any function A at temperature T is given by 

= ^W, (2) 

where Tr is a shorthand for the sum over all the spins = ±1 , and (3 = l/k^T is the 
inverse temperature. 

We define the correlation length in the horizontal direction £i(T) by 

£i(T) = lim —rr. (3) 

— m\CT00^2™0/T 

Similarly, we define the correlation lengths £2 CO and ^(T) for the other two directions as 
in Fig. EC 

Our definition ([3]) of the correlation length is motivated by the observation that the 
correlation between spins at two ends of the lattice should behave as (o" ocr 2 n )^ , ' +1 ' > pa 



3 



000 ~ a/o 
* 2" * 

Figure 3: £i(T), £2 CO and ^(T) are the correlation lengths in the three directions. 



exp[— 2 n /^i(T)}. This is apparently different from the usual definition, motivated by the 
observation that (<J x <J y )T ~ exp[— \x — y\/£] for any sites x, y. For the Ising model on a reg- 
ular lattice, however, the two definitions are known to give the same correlation length (at 
least in the high-temperature region). Here we conjecture equivalence, and therefor employ 
the simpler definition, (j3J). 

We can now state our main results. We choose the interactions Ji, J 2 and J 3 so as to 
satisfy < J\ < J 2 < J 3 . When the temperature approaches zero (which is the critical point 
of the model), the correlation lengths diverge. For the ratios of the correlation lengths, we 
show that as T — > we have 



6C0 6C0 
a(T) ~" ' ^(r) 



oo (4) 



if 



£>^ + l, (5) 



which corresponds to the shaded region 1 in Fig. HJ and 

feCO 6(r) 

' 6 CO 

if 



(6) 



< -p + 1, (7) 
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Figure 4: The region 1 is where J3/J1 > (2J 2 )/J\ + 1, and region 2 is where J3/J1 < 
(2J 2 )/Ji + 1. We have partial and complete restoration of anisotropy in these two regions, 
respectively. 



which corresponds to the shaded region 2 in Fig. |H These correspond to the partial and 
perfect restoration of isotropy discussed in the introduction. 

Moreover, when the interactions satisfy (j3J), and we have a partial isotropy restoration 
(jl]), we can show that the leading-order singular behavior of the correlation lengths is given 
by 



f 2 (T) ~exp (ln2)e 



3 2(Ji+J 2 )/0 B T) 



cxp 



(ln2)e : 



2(Ji+J 2 )/(k B T) 



cxp 



2(1/3 — Jl ~~ 2^2 



(9) 



as T — > 0. We see that the correlation lengths diverge as a double exponential of l/(k B T), 
while the ratio ^ 3 (T)/^i(T) diverges as a single exponential of 1/(/cbT). (The double expo- 
nential behavior is peculiar to the Ising model on the SG, and it was previously found in the 
isotropic model with J\ = J 2 = J3. See Ref. [8].) 

In the next section, we derive these results using the exact renormalization group method. 
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Figure 5: Thanks to the special geometry of the SG, the renormalization group map can be 
obtained if we consider a decimation procedure for the second generation SG. 

3 Derivation 

3.1 Renormalization group transformation 

It is well known that various models on the SG can be treated using the exact renormalization 
group transformation (usually called "decimation"), in which one sums up parts of the 
degrees of freedom. 

Let 7i n be the Hamiltonian of the form ([I]) on the n-th generation SG, and let ^(a^)] 
be an arbitrary function of a xy with both x and y even. We denote by Tr the sum over 
a xy = ±1 for all the sites xy with both x and y even. Let Tr" denote the sum over a xy = ±1 
for the rest of the sites. Then from (j2J), we find that 

{A) (p) = LlAe _ = IrAlr e_ = lr Ae ^ 
Tre~ n " Tr Tr e Wn Tr e _w « 

where we have set H n = (37i n and defined H' n by exp (— Ti'A = Tr"exp (—Unj. Note that 
the final expression can be regarded as the thermal expectation value in the Ising model on 
the (n — l)-th generation SG with the Hamiltonian H! n . 

Because of the special geometry of the SG, the interactions in the new Hamiltonian (as 
functions of the interactions in Ti,' n ) are independent of the generation n. Therefore, it suffices 
to consider the case with n — 2. 

Let H be a Hamiltonian on the second generation SG, and write it as 

H = f3H = -Ki(cr 1 a 2 + a 2 & 3 + 030-1) 
-K 2 (cr 1 a 2 + a 2 a 3 + a^) 

-K z {a x o 2 + a 2 cr 3 + a 3 ai), (11) 

where Ki = /3Ji,i = 1,2,3 (see Fig. [5]). (Below, we use Ki instead of /3Jj.) Because we 
have Ki oc 1/T, T — > implies Ki — > 00. The renormalized parameters K[, K' 2 and K' 3 are 
determined from 

E E E e-e» = e-e* (12) 

cri=±l fT2=±l CT3=±1 
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and 

V! = -K' x o x (y 2 - K' 2 cr 2 cr 3 - K3O-3O-1. 
An explicit calculation shows that 

Ecosh^i) 
Acosh(2K 2 )cosh(2i\:3) 

J Bcosh(2iT 2 ) 
Acosh(2 J fsT 3 )cosh(2ii'i) 
Bcosh(2K 3 ) 











K' 2 = 


ii„ 
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K = 


iln 
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v4cosh(2iri)cosh(2fs: 2 ) 



(13) 

(14) 
(15) 
(16) 



where 



and 



A = exv^ + Kz + Kd+^exvQKi-Ki-Ki-Ks), 

8=1 



B = exp( J ftT 1 + K 2 + AT 3 ) coshp^ + K 2 + K 3 )\ 
3 

+ £ exp{2Ki -Kx-Ki- K 3 ) cosh [2 (^ + K 2 + K 3 - 2Ki)\. 

i=l 



(17) 



(18) 



We set Ki(0) = (3Ji for % = 1, 2, 3. We apply the renormalization group (RG) map given 
by Eqs. (ITlft - flEl) to the initial conditions Ki(0), K 2 (0) and -^3(0) repeatedly. We denote 
the results of n applications of the RG map by Ki(n), K 2 {n) and K 3 (n). 



3.2 Correlation lengths 

Here we derive a general expression for the correlation lengths using the RG method. Let 
7i n be the Hamiltonian of the form (JT]) on the n-th generation SG. Then, using Eq. ffTQ]) n 
times, we find 

(o"ooc r 2"o)f +1) = (o"ooO"io)5Jj, (19) 

where the right-hand side represents the expectation value in the model on the first genera- 
tion SG, and 

H[ = -K 1 (n)a 00 a 01 - K 2 {n)a m a w - K 3 (n)a w cr 00 . (20) 
By explicitly calculating the right-hand side of Eq. ffl9|) . we find 

, (1) _ e^'W cosh (K 2 (n) + K 3 (n)) - e~ K ^ cosh (K 2 (n) - K 3 (n)) 

{a ° oCTlo) ni ~ e K,(n) cosh (^ 2 ( n ) + ^ 3 ( n )) + e - Kl {n) cosh [x^n) - K 3 \n)) ' { ' 



To evaluate the correlation lengths, we let hq be a sufficiently large number such that 
Ki(n) <C 1 holds for all n > no- Then, for n > n , the RG map given by Eqs. (fl4|) - 
( JT6l) is approximated by 

+ 1) ~ ^(n) 2 , (22) 

whose solution is 

K,(n)^^K) 2n ""°, (23) 

for % = 1, 2, 3. 

When quantities Ki(n) are all sufficiently small, the right-hand side of Eq. (121]) is ap- 
proximated by Ki(n). Then, from the definition ([3]), we immediately get 

Exactly the same analysis is used for the other directions, and we get 

2 n ° 

« r » " -i^K) < 25 > 

for 2 = 1,2 and 3. 



3.3 Behavior of the correlation length in the isotropic case 

Before dealing with the problem of isotropy restoration, we briefly examine the singular 
behavior of the correlation length in the isotropic model with J\ — J2 — Jz- (Because Ki(n) 
are identical for i — 1, 2, 3, we drop the subscript i hereafter.) Since we have the expression 
(jSJ), it suffices to evaluate no- 

As we are interested in the T — > limit, we assume K(0) 3> 1 and note that the quantities 
in Eqs. f|Til) - ffTBl) are approximated as 

K(n + K{n) - exp(-4K(n)) (26) 

if K(n) ^> 1 holds. This difference relation is approximated accurately by the differential 
equation 

4-K(n) = -e 4i ^\ (27) 
an 

whose solution is 

K{n) = ~ In (e 4iY(0) - An) . (28) 

To check the accuracy of the approximation we compare in Fig. |6] the RG flow obtained by 
numerically solving Eqs. (|T^I) - (fT6l) with the solution (1281) . We find that the approximation 
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Figure 6: Each diamond denotes 500 points of the renormalization flow, while the curve is 
the approximate RG flow of Eq. (I28I) . 

([28]) is quite accurate as long as K(n) > 1. Furthermore, we observe that K(n) immediately 
enters the region K(n) <C 1 when the condition K(n) ^ 1 is violated. Based on the second 
observation, we assume that the behavior of the RG flow in the regions K(n) <C 1 and 
K(n) ^> 1 can be connected almost smoothly. 

To determine n , we set K(n ) = K , where K <C 1 is an arbitrary constant that we 
shall fix. Then, from Eq. (I28j) . we get 



,4if(0) 



3 4Ko 



By substituting this expression into Eq. (1251) . we finally get 

2 I( e 4X(0)_ e 4ir ) 



(29) 



e(r) 



cx exp 



const, e 



4K(0) 



cxp 



Ce 4.J/(k B T) 



(30) 



in agreement with Ref. [8]. Here, C is a constant (which cannot be determined precisely 
from the present crude analysis). 

3.4 Partial and perfect restoration of isotropy 

We now treat the general case with Ki(0) < K2 (0) < -^3(0). We concentrate on the region 
in which Ki(n) ^> 1. The RG flow in this region essentially determines the low temperature 
behavior of the system. Here, the RG map (|l4p is approximated by 



4 



e ~2(K t (n)+K 2 (n)) + e -2(K 2 (n)+K 3 (n)) 



-2(Jf 3 (n)+ifi(n)) 



e -4Ki(n) + e -4K 2 (n) _|_ e -4K 3 (n) 



(31) 
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Figure 7: The RG flow in the two-dimensional space (K^/K^^K^/K^) obtained from 



(EH). There are three fixed points, (K^ n) /K[ n \ K^/K[ 



l,oo), (1,3) and (1,1) 



We get similar equations for K 2 and K 3 from Eqs. (fl5l) and ({16]) . respectively. (They are 
obtained by cyclic permutations of the indices 1, 2 and 3.) 

Let us first focus on the behavior of the RG flow on the line K 2 (n) / Ki(n) = 1 and further 
assume that K^{n) — Ki(n) ^> 1. Then the RG map fl3T|) is approximated by 



K*(n + 1) 



K 3 {n) 

KAn + 1) K x (n) 



kK x {n)e 



-1 KM 



Ki(n) 



(32) 



This implies that in the two-dimensional space, the point (K2(n) / Ki(n) , K 3 (n) / Ki(n)) = 
(1, 3) is a fixed point. Moreover, we find that there are RG flows which emerge from the 
fixed point, as depicted in Fig. [7J 

Next we consider the case with K 2 (n) / Ki(n) > 1, which is the region 1 in Fig. |H From 
Eq. (J3ID we find that K 2 (n + 1)/Jfi(n + 1) > K 2 (n)/K 1 (n) and the RG map (EH) is 
approximated by 



K 3 (n) 
KM 



+ 1 



' KM 

km 



+ 1 



' K 3 (n + 1 
KUn + 1 




K 2 {n + 1) 
#i(n+l) 



+ 1 



(33) 



It is hypothesized that the RG flow in the case with K 2 (n) / Ki(n) > 1 is like that depicted 
in Fig. [71 In the case with K 3 (n)/Ki(n) > 2K 2 [n) / Ki[n) + 1, the RG flow goes to the 
fixed point (K 2 (n) / Ki(n) , K 3 (n)/Ki(n)) = (l,oo), and in the case with K 3 (n)/Ki(n) < 
2K 2 (n)/K 1 (n) + l,th.e RG flow goes to the fixed point {K 2 {n)/ K i{ n )i K s (n)/KM) = (M)- 
We confirm this picture in each case in what follows. 
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3.4.1 Partial restoration of isotropy 

Now we treat the case with K%(n)/K\(n) > {2K 2 {n)) / K\{n) + 1. We separately consider 
three different regimes of n. Let us define n\ by the relations K\{rii) ~ K 2 {nx) — 1, -^3(^1) ^> 
1. (Below we confirm that such an n\ exists.) Then Eq. ( 13T1) transforms into 

dK l( n ) _ _I ^ e ~2(X 1 (n)+^ 2 (n)) _ e -4Ki(n) + e "4X 2 (n)^ ? ^ 



(in 4 

„ _I |' e -2(^ 1 (n)+^ 2 (n)) + e -4Jfi(n) _ g-4X 2 (n) 

(in 4 ^ 



(35) 



„ _1 (' e -2(^ 1 (n)+^ 2 (n)) + e -4Jfi(n) + g-4X 2 (n)\ /ggN 

(in 4 V / ' v ; 

where we have used the relation K 3 (n) — Ki(n) ^> 1. These equations can be solved analyti- 
cally Letting L x [n) = e-^W+^W), L 2 (n) = e ~^K 2 (n)+K 3 (n)) and L ^ = e -2(x 3 (n)+K 1 (n)) ) 

the equations (I34"I) . (135]) and (!36|) become 

~ Lin , 37 

n 

dL 2 (n) 



n 

and 



Li(n)(L 2 (n) + L 3 (n)) (38) 



respectively. Then, from Eq. ( |371) . we have 



=L 1 (n)(L 2 ( n ) + L 3 ( n )), (39) 



- t^§k (40) 



Next, we define L±(n) = L 2 (n) ± L 3 (n), and we have 



d£+(n) r , w . , 2Li(0) . . 

2Li(n)L+(n) = v / L+(n). (41) 



iV ' ^ v ' 1-Li(0)n 
Solving this differential equation, we get 



Mf0 * ^a+^a. ,42) 

+ V ; (l-Li(O)n) 2 V ; 



We can obtain L_(n) similarly, using 



dL_{n) 
dn 
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(43) 



from which we have 



From Eqs. (1421) and 

L 2 (n 



L 3 (n) 



L_(n)~L_(0) = L 2 (0)-L 3 (0). 
we find 

L+(n) + 
2 

L 2 (0) + L 3 (0) + (L 2 (0) - L 3 (0))(1 - Lx(O)n) 



(44) 



2(1 -Li(O)n) 2 



L_i_(n) — L_fn) 



L 2 (0) + L 3 (0) - (L 2 (0) - L 3 (0))(1 - M0)n) 
2(1-Li(0)n) 2 

Therefore Xi(n), K 2 (n) and X 3 (3) are finally expressed as 

1 



In 

4 

I In 

4 



In 



£ 2 (n) 


Li 


(n)L 3 


[n)_ 




Ls(n) 




Lx 


(n)L 2 


[n)_ 




Li(n) 




L 2 


(n)L 3 


[n) 



From Eq. (1471) . we see that ni can be taken as 

m ~ e 2 ( L i(°)+ L 2(°)). 



(45) 

(46) 

(47) 
(48) 
(49) 

(50) 



We can confirm that i^i(ni) ~ K 2 (n 1 ) ~ 1 by substituting n x into Eqs. (j4"Tj) and (j4"8~l) . 

In Fig. El we compare the approximate RG flows (j4Tj) -( l4"9l) with numerical data. The 
agreement is again quite satisfactory. 

Next we consider the regime in which n 2 > n > ni, in the case that n 2 is such that 
K x (n 2 ) < l,fG(n 2 ) < 1 and K 3 (n 2 ) ~ 1 . In this regime Eqs. (jI5j) and (JTSD are 

approximated as 



+ ~ 4Ki (n), 
K 2 (n+1) ~ 4K 2 (n), 

# 3 (n+l) ~ K 3 (n)-^ln2. 



The relations (1511) . (152]) and (153]) are solved as 

i^(n) ~ (2 n ~ m - l)i^ 1 (n 1 ) 2 " - ' 11 . 



(51) 
(52) 

(53) 
(54) 
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Figure 8: Each square, plus sign, and diamond represents 2000 points of K±, K 2 and K 3 , 
respectively. The curves are the corresponding approximate RG flows of Eqs. (]4TT ) — ( ]49]) . 

K 2 (n) ~ (2™-l)K 2 (n 1 ) 2 """ 1 , (55) 
K 3 (n) ~ K 3 (m)- (n ~ ni) ln2. (56) 

From Eq. fl56l and K 3 {n 2 ) — 1, n 2 is estimated as 

n 2 ~ m + ^- (K 3 (0) - K^O) - 2K 2 (0)) . (57) 
In 2 

The correlation lengths are estimated by using Eqs. (12"5"|) and (I5"l"j) - (l5"6"j) : 

W) ~ -ln^K) (58j 
^ 2 cxp[2(K 1 (o)+x 2 (o))] ) (59) 

£ 2 (y) ~ 1 2 ex P' 2 < K i(°>+^(o))l (60) 

-lni^ 2 (ni) 

^ 2 cxp[2(^ 1 (0)+A' 2 (0))] ; (g^ 
~ 1 2Cxp \ 2 (K 1 (0)+K 2 (0))] + T ^(K s (0)-K 1 (0)-2K 2 (0)) 

-In If 3 (no) 

_ 2exp[2(ifi(0)+K 2 (0))]+ I | 5 (if3(0)-K 1 (0)-2Jf 2 (0))_ ^ 

The ratio of ^(T) to £i(T) [or ^(T 1 )] diverges as T — > oo. 
3.4.2 Perfect restoration of isotropy 

Finally, we show that in the region 2, we have K 2 /K\ — > 1 and K 3 /K\ — > 1 as T — > 0; i.e, 
there is perfect restoration of isotropy. We treat the following three cases separately. 
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Figure 9: The region 2a with K 2 /Ki > 1, K 3 /K 1 < 3, and the region 2b with K 3 /K 1 < 
2{K 2 /K l ) + 1, Ks/K, > 3. 




Figure 10: The RG flow in the two-dimensional space (fTi(n), K 3 (n) / Ki(n)) . 
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(i) The case K2/K1 = 1. 

Assuming K 3 (n) — Ki(n) ^> 1 in the RG map (I3TT) . we get 



K 3 (n) K 3 (n + l) _ e 4 * 1 ^ f AT 3 (ti 



#i(n) #1(71+1) AK^n) { Kx{n) 

and 

e -4iCi(n) 



(64) 



K x {n)-K x (n + \)~ . (65) 



Then, with x n = Ki(n) and y n = Kz(n)/K\(n), these two relations yield 

2/n+l — Un 3 — ?/ n 



(66) 



This becomes the differential equation dy/dx = (3 — y)/x whose solution is y = 3 — C/x, 
with an arbitrary constant C. Therefore, K\(n) and Kz(n)/K\(n) vary on the curve 

as in Fig. [TUJ When n is large and the assumption K 3 (n) — K\ (n) ^> 1 is no longer valid, the 
RG flow deviates from the curve (|67|) . Let 1 + e be the value of K 3 (n) / Ki(n) at the value 
of n where this deviation starts to take place. We can take e as small as we wish by letting 
T — > (i.e., Ki(0) — > 00). Therefore, we conclude that K3/K1 — > 1 as T — > 0. 

(ii) The cases K2/K1 > 1 and K3/K1 < 3 (the region 2a in Fig. 

As in the case (i), the RG map (|3"T|) and the assumption K${n) — K\(n) ^> 1 yield the 
equation (I64"l) . The equation (|65|) is modified as 

-4Xi (n) 

^(n) -K^n + 1) < , (68) 

which is a stronger relation than Eq. fl65l) . Therefore, the same argument as in (i) shows 
that K 3 (n) / Ki(n) can be made arbitrarily close to 1 by letting T — > 0. 

fzzzj TTie case #3/^1 < 2(K 2 /K 1 ) + 1 and # 3 /#i > 3 (Ue region 2b in Fig. \®. 

From Eq. (j32"I) . we find that the RG flow depicted in the {Ki/K\, K3/K1) plane follows 
the curves in Fig. [9j Thus the RG trajectories flow from the region 2b into the region 2a. 
In what follows, we show that the conditions K^{n) — Ki(n) ^> 1 and K2(n) — K\(n) ^> 1 
are satisfied when the trajectory enters the region 2a, provided that we take a sufficiently 
small values of T. Then we can use the argument given in (ii) to conclude that a perfect 
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restoration of isotropy takes place for T — > 0. Supposing that K 3 (n) — K\(n) 1 and 
K 2 {n) -K x {n)^>l, the RG map yields 

p -4Ki(n) 

^(^-^(n+l) ~ . (69) 

Because this shows that K\{n) increases with n, we see that the desired condition K 3 {n) — 
Ki(n) ^> 1 is valid when the trajectory enters the region 2a. Let a be the value of 
K 2 (n) / Ki(n) at which the trajectory enters 2a. By taking Ki(0) so that (a — 1)X 1 (0) ^> 1, 
we can guarantee that the desired condition /^(n) —K\{n) ^> 1 is valid when the trajectory 
enters 2a. 

3.5 The differences between the present and previous studies 

Brody and Ritz [9] , and later Bhattacharyya and Chakrabarti [TD] studied the same problem 
of the restoration of isotropy in the Ising model on the Sierpihski gasket. They both started 
from a partially anisotropic model with J 2 = ^3 or J\ = J 2 , and studied the model using 
the same RG equation. But the conclusion of Ref. [9J is that there is complete isotropy 
restoration while that of Ref. [10J is that isotropy restoration is impossible. Note that our 
conclusion that the restoration may be either partial or perfect agrees with neither of these 
two conclusions. We stress that our analysis is much more elaborate and careful than the 
previous two, and it therefore provides the most reliable conclusion. 

Let us point out some of the inaccurate points in the previous works, which may be the 
sources of (partially) incorrect conclusions. 

In Ref. [10J, the authors make a rather crude "weak coupling" approximation to the 
RG map, which effectively brings the system to an extremely high temperature by a small 
number of application of the RG map. Thus is is not possible for nonlinear effects to exist, 
and hence there results no isotropy restoration. 

In Ref. [9], the authors make an unjustified approximation using a philosophy very 
different from that in Ref. [10J. When deriving their key equation, (7), the authors rely on a 
"rescaling" argument which can never be exact. Consequently, (7) becomes an approximation 
which is probably valid only when K n+ i ~ 1. But this condition cannot be assumed for 
L n < 1. 

In fact, our analysis clearly demonstrates that the present problem cannot be under- 
stood using a single approximate RG map. Instead, one needs to carefully use different 
approximations depending on the ranges of the parameters, as we have done here. 
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